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ABSTRACT 

Inspection of broadband VLF records from the Stanford University/ 
Stanford Research Institute experiment aboard 0G0 4 shows that some 
whistlers exhibit a banded structure in which one or more bands of fre- 
quencies are missing from the whistler's spectrum. The phenomenon is 
commonly observed by satellites on midlatitude field lines at all local 

times and at various longitudes around the world. The dispersion of 

1/2 

banded whistlers (BW) is of several tens of sec , indicating that they 
originated in the opposite hemisphere and are propagating downwards at 
the satellite. BW are generally spread in time (tenths of seconds) 
rather than sparply defined and tend to occur at random. The frequency 
spacing of the bands may be either uniform or irregular, and may vary 
radically between successive events. Several possible explanations for 
BW are considered. In particular, an analysis of the interaction of plane 
electromagnetic waves traveling in an anisotropic plasma with a field 
aligned slab of enhanced ionization is presented with promising results. 
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I . INTRODUCTION 


A. INTRODUCTION 

A detailed study of broadband very low frequency (VLF) whistler 
spectra reveals fine structure in virtually every case. Perhaps the 
most spectacular type of fine structure is the phenomenon of "banding” 
in which one or more bands of frequencies are missing from the whistler's 
spectrogram. A particularly striking example is shown in Figure 1. The 
objective of this report is to describe and explain observations of this 
phenomenon . 

B . BACKGROUND 

The surge of current which occurs during a lightning flash radiates 
a spectrum of electromagnetic energy which extends from a few Hertz up to 
around 20 MHz. The spectrum of these radiated waves generally peaks in 
the range 5-10 kHz. These waves may eventually enter the ionosphere and 
magnetosphere. At VLF frequencies the properties of the medium are such 
that the waves follow curved paths which are roughly aligned with the 
earth's magnetic field. Depending on the angle which the wave normal 
direction makes with the magnetic field the wave may also move in or out, 
to lower or higher L shells as it travels. At times such waves may become 
ducted (trapped by small field -aligned irregularities in the plasma density) 
and thus will more faithfully follow a particular field line. In addition 
the medium is dispersive so that different frequency components travel at 
different velocities. In the audio range this often results in the 
production at a receiver of a descending tone, called a "whistler.” 

Whistlers were first explained by Storey [1953]. Since then many 
publications on whistlers have appeared. For more details on the phenomenon 
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OGO 4 (ACT) 
03 OCT 67 


f 

kHz 

10 - 




FIGURE 1. FREQUENCY-TIME SPECTROGRAM OF A BANDED WHISTLER (BW) SHOWING 

AN UNUSUALLY UNIFORM BANDED STRUCTURE. As usual, the dispersion 
of this event (of the order of 100 sec 1 / 2 ) indicates that it originated in 
the opposite hemisphere and is downgoing at the satellite. The relative 
darkness of the display measures the relative strengths of the received 
signals. The nulls in the spectrum of this banded whistler occur at odd 
half multiples of 810 Hz. The ephemeris indicates a proton gyrofrequency 
of 719 Hz and an invariant latitude A - 51° at the satellite. The local 
time is a few minutes after midnight. Also visible at the left of the 
record is an SP whistler which displays similar, but not identical banding. 
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see Helliwell [1965]. 


The earliest published examples of banded whistlers known to the 
author were from ground whistler data in Helliwell [1965]. However, 
examples of banded whistlers in the ground data are extremely rare. The 
present study concentrates on data obtained with the Stanford University/ 
Stanford Research Institute VLF experiment flown on 0G0 4, a low altitude 
polar orbiting satellite. 

The Orbiting Geophysical Observatory, 0G0 4, was launched on 28 July 
1967 into a polar orbit with an inclination of 86°. The orbital period 
is around 100 minutes and the height varies from approximately 400 to 900 
km. The satellite carries more than 20 different experiments. The data 
for this report were obtained from experiment D02, the Stanford University/ 
Stanford Research Institute VLF experiment . 

Experiment D02 consists of a 2.9 meter diameter loop antenna and a 
broadband preamplifier that drives three narrow band stepping receivers, 
two broadband receivers and one phase tracking receiver. The antenna was 
designed to be switchable on command so as to act as either an electric 
or magnetic field sensor. However, most of the data were obtained with 
the magnetic sensor because of a malfunction in the experiment. Table 1 
describes the characteristics of the two receivers which provided data 
for this report. For a more detailed description of the experiment see 
Rorden et al [1966]. 

In Chapter 2 data from the above experiment which characterize the 
banded whistler phenomenon are presented. Chapter 3 will present some 
possible explanations, the most promising of which, spectral reflection 
or transmission from a field aligned slab of enhanced ionization, will 
be examined in detail. A summary, discussion, and suggestions for further 

research are presented in the final chapter. 
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TABLE 1. RECEIVER CHARACTERISTICS OF EXPERIMENT D02 



Narrow bands 

Broadbands | 

Band 1 

Band 2 

Band 3 

Ph .Tracking 

VLF 

ELF 

| Frequency range (kHz) 


BBH 

12.5 - 100 

14.4 - 26.3 

0.3-12 .5 

0.015-0.300 

Threshold 

Sensitivity 

magnetic field (piy) 

32-4 

6 - 0.8 

1 . 2 

1.6 

8000-8 

10 5 - io 3 

electric potential 
at input of pre- 
amplifier (fiV) 

0*6 - 0 *06 

0.06 

0.08 

0.06 

80 

200.0 - 20 

Dynamic range (db) 

90 

90 

90 

80 

80 

90 

Local oscillator freq<> (kHz) 

3.16 - 4.53 

25.33 - 36.25 

202.5 - 290 

23 - 34 and 
2.7 - 3.6 



IF freq. (kHz) 

2.97 

23.75 

190 

8 . 2 and 5 . 0 



3-db bandwidth (Hz) 

40 

150 

500 

50 



AF/step ( Hz/ step) 

5 .4 

43 

344 

100 


l 

Data rate (kbit/sec) 

4,16,64 

4,16,64 

4,16,64 

4, 16,64 



Sweep rate (steps/sec) 

3.36,13.9,56 

3.36,13.9,56 

3.36,13.9,56 




Sweep time (sec/sweep) 

73.4,18.4,4.6 

73.4,18.4,4.6 

73.4,18.4,4.6 




Integration time constant (msec) 

460,114,20 

146,40,9.4 

48.4,18,4.4 

32,8.6,2 

























































II. DATA 


Broadband VLF spectra show a variety of fine structure, some of 
which are illustrated in Figure 2. Figure 2a shows a sharp well defined 
whistler which has gaps in its frequency spectrum. Close inspection of 
the data reveals such gapping in most well defined whistlers. A somewhat 
rarer phenomenon is illustrated in Figure 2b which shows a whistler having 
peaks of enhanced amplitude in its spectrum. This "beading" effect is no 
better understood than the spectrum gaps of Figure 2a. The fine structure 
of Figures 2a and 2b may be related to the banded whistlers. A third type 
of fine structure is illustrated in Figure 2c, where it is seen that the 
trace is composed of a collection of "blobs" rather than a single trace. 
These "blobs" are generally of a size much greater than the resolving limits 
of the spectrum analyzer. In this case the entity we call a whistler 
appears to consist of a collection of separate and distinct wave packets 
each with its own set of frequencies and propagation paths. 

Banded whistlers (Figure 1 is a particularly outstanding example) 
are further illustrated in Figure 3, which shows six minutes of data taken 
over Orroral, Australia (ACT) on 30 May 1968. While Figure 1 shows a highly 
uniform spacing, this pass shows quite different behavior. For the whistler 
labeled "e", for example, the spacing between nulls ranges from 430 Hz to 
1000 Hz. Also, little correlation is found from one banded whistler to 
the next. In "b" the spacing ranges from 300 to 500 Hz while 18 seconds 
later ("C") a relatively uniform spacing of approximately 1300 Hz occurs. 
Proton gyrofrequencies, deduced from the ephemeris, range smoothly from 
693 Hz at "A" to 728 Hz at "b" . 

The phenomenon of banding occurs worldwide, having been detected 

at all receiving stations. Out of a sample of 87 passes one or more banded 
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OGO 4 JOB 16 DEC 67 -1600 UT 
2a 2b 





FIGURE 2 * FREQUENCY -TIME SPECTROGRAMS ILLUSTRATING THREE TYPES OF 
WHISTLER FINE STRUCTURE. a) Shows a sharp, well defined 
whistler with two of the gaps in its spectrogram indicated by arrows, 
b) Shows a whistler having peaks of enhanced amplitudes, c) Shows a 
broadly dispersed whistler trace. The arrows indicate some of the small 
discrete elements of the whistler trace. 
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OGO 4 (ACT) 
30 MAY 1968 


f 

kHz A - 35°, f b = 640 Hz 



1410 UT 

A -39°, f b -669 Hz 



14)1 

A = 43°, f b =693 Hz 



1414 

A = 55°, f b =745Hz 



1415 


FIGURE 3. FREQUENCY -TIME SPECTROGRAM OF BROADBAND DATA FROM OGO 4 SHOWING MANY EXAMPLES OF BANDED 
WHISTLERS. Some of the more prominent examples are labeled "A" through "e" . 


whistlers occur on 70% of the nighttime passes. It occurs on 35% of the 
daytime passes, but this difference appears to be wholly attributable to 
the reduced overall whistler activity in daylight hours (i.e., on daytime 
passes there is at least an order of magnitude fewer whistlers than on 
nighttime passes — in fact often less than a dozen whistlers are detected 
per pass) . Banded whistlers are well represented in those daytime passes 
which detect reasonably high levels of whistler activity. 

At this point it is reasonable to ask "What is the distinction 
between the phenomena shown in Figures 2a and 2b and the banded whistlers?” 
At Stanford University our standard method of data presentation involves 
a 0-10 kHz or 0-13 kHz spectrogram displayed on a 3-cm wide roll of film. 

On this scale the phenomena shown in Figures 2a and b tend to be observ- 
able only on close inspection. By virtue of the number of nulls and/or 
the increased frequency width of the nulls, that subset of the data referred 
to as 'banded whistlers” tend to be strikingly obvious to a casual observer. 
While this is admittedly a subjective distinction, it is operationally 
effective since the number of borderline cases is too small to signifi- 
cantly affect the statistics in the preceding and following paragraphs. 

In addition to OGO 4 the phenomenon of banding has also been observed 
in this laboratory's data from OGO T s 1, 2, and 3. Because of the availa- 
bility of large quantities of data all statistics are based on 0G0-4 data. 
Because of the high rate of banded whistler activity at ACT and the fact 
that ACT 'sees” poleward beyond the plasmapause and to very low latitudes 
in the equatorial direction, more effort was spent analyzing data from 
this station than any other. 

Figures 4 and 5 show the occurrence of banded whistler activity as 
a function of satellite dipole latitude and L value respectively. The high 
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GEOMAGNETIC LATITUDE 


FIGURE 4. HISTOGRAM SHOWING THE OCCURRENCE OF BANDED WHISTLERS AS A FUNCTION OF DIPOLE LATITUDE 

(JENSEN-CAIN MODEL) OVER ORRORAL, AUSTRALIA. The data cover the periods September 1-3, 
18-19, 1967, October 2-4, 1967 and May 26-30, 1968. These periods were selected because they each 
contain at least one pass with many well defined banded whistlers. Approximately 80% of these 
whistlers occurred at night. Data were generally received from the satellite beyond the upper and 
lower latitude cutoffs in the occurrence of banded whistlers. 


o 



NUMBER OF 
BANDED WHISTLERS 


W 

t - 1 


I 

0 

01 



FIGURE 5. SAME DATA AS PRECEDING FIGURE DISPLAYED AS A FUNCTION OF L-VALUE 


latitude cutoff in these figures appears to coincide with the crossing of 
the plasmapause since banded whistler activity has been observed right up 
to the point where all whistler activity disappears. The low latitude 
cutoff at ~ 30° geomagnetic latitude appears to be real, although the 
sharpness of the cutoff in Figure 5 is rather misleading since, due to 
geometrical factors, the satellite spends greater amounts of time at the 
lower L shells. Data from the other stations are consistent with these 
latitude variations . 

From records such as Figure 1 it is clear that severe attenuation 
of the signal has occurred in the nulls. The broadband receivers used on 
OGO 4 obtain their large ( ~ 60db) dynamic range by log compressing and 

clipping the received signals. In order to accurately determine the amount 
of attenuation in the nulls it is necessary to use the data from the 
narrowband stepping receiver (calibrations are shown in Figures 6 and 7) . 

The band-2 narrowband receiver sweeps upward in frequency at 43 Hz/step 
at a specified number of steps/second depending on groundbased command 
signals. It is rare for a well-defined null at some frequency in a banded 
whistler 'to occur just as the stepping receiver is scanning across that 
particular frequency. However, two clear examples have been found where 
the stepping receiver sweeps in frequency from a maximum to a null in a 
banded whistler. Figures 8 and 9 show the sweeping receiver data along 
with the corresponding banded whistlers. With the aid of the calibration 
charts Figure 10 was plotted showing the observed signal strength recorded 
by the sweeping receiver as the sweep moves into a null in the corresponding 
banded whistler. In both whistlers attenuations of — 14 db were noted. 
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OUTPUT (volts) 



FIGURE 6- OGO-4 CALIBRATION OF BAND 2 RECEIVER IN MAGNETIC 
MODE AT 10 kHz, Use Figure 7 for calibration at 
other frequencies. 
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DIGITAL OUTPUT 




























OGO 4 ACT 3 OCT 67 


OUTPUT 

VOLTS 


f 

kHz 


5 - 


BAND I 
O- 
5 - 



BAND 2 


O- 


BAND 3 


| 1452: 40 UT 



I SECOND 


FIGURE 8 - COMPARISON OF QGO-4 BAND 2 AND BROADBAND RECEIVER DATA AT 
- 1452:40 UT ON 3 OCTOBER * 1967 OVER ORRORAL, AUSTRALIA. 

This record includes the data shown in Figure 1. The top graph shows 
the amplitude vs time variation of the signal received by the sweeping 
receivers. The Band 2 receiver has a 150 Hz (3 db) bandwidth and sweeps 
up in frequency at 43 Hz/step, 13.9 steps/second, repeating the sweep 
every 18.4 seconds. The bottom panel shows the frequency vs time spectrum 
of the data from the broadband receiver. The upward sloping straight line 
in this figure indicates the midfrequency being observed at any particular 
time by the Band 2 receiver. The time scales are the same on both panels. 
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OGO 4 ACT 3 OCT 67 


OUTPUT 

VOLTS 


f 

kHz 


5 - 


BAND 1 


o- 

5- 



BAND 2 


°1 


BAND 3 


| 1315:30 UT 



I SECOND 


FIGURE 9. 


SIMILAR TO FIGURE 8 EXCEPT THAT THE OBSERVATIONS 
WERE MADE ON THE PRECEDING PASS AT ~ 1315:30 UT . 
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FREQUENCY, kHz 


FIGURE 10. DATA FROM FIGURE 8 CALIBRATED FROM 

FIGURES 6 and 7 TO SHOW THE SPECTRAL 
AMPLITUDE VARIATION IN THE VICINITY 
OF THE 2.5 kHz BAND IN FIGURE 1. 
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III. 


INTERPRETATION 


A number of potential explanations for the banded whistlers can be 
advanced. The merits of eight models are discussed in this chapter. 

They fall roughly into the following five categories: 1) instrumental 

effects; 2) superposition of multiple signals; 3) banding in the source 
spectrum; 4) gyroresonance interactions; and 5) specular reflection or 
transmiss ion . 

A. ANALYSIS OF POTENTIAL EXPLANATIONS 
- 1 . Instrumental Effects 

It is, of course, necessary first to deal with the possibility 
that the phenomenon is spurious. Potential sources of instrumental 
effects are a) the satellite; b) the receiving station, c) the processing 
at Goddard Space Flight Center, and d) our Rayspan spectrum analyzer. 

The occurrence of the phenomenon on other satellites (e.g., 

OGO 2, T. Bell, private communication) and all ground receiving stations 
argues against (a) or (b) . The apparent latitude effects noted in 
Figures 4 and 5 argue very strongly against (b), (c), or (d) and to a 
lesser extent, (a). The possibility that non-uniform response of the 
Rayspan filters is responsible for the phenomenon is unreasonable in view 
of the randomness of the frequencies at which the nulls occur. This source 
was further checked by comparing the same data on 0-10 kHz and 0-13 kHz 
records. The Rayspan (spurious) lines remained fixed while the nulls 
moved to their new positions, as expected. Finally, ground whistler data 
processed on the same Rayspan analyzer have only rarely contained any 
banded whistlers. 

2 . Superposition of Multiple Signals 
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Multiple stroke. Given that the effect is real, one possibility 


is that the nulls are caused by destructive interference of signals from 
multiple stroke lightning. Assuming lightning to be impulsive the total 
signal V(t) is given by 


V(t) = V 1 6(t-t 1 ) + V 2 6(t-t 2 ) 


( 1 ) 


The Fourier spectrum of V(t) is given by 
V(tD) = 


V(t) dt = V e ~ iCDt l + V Q e" iCUt 2 


( 2 ) 


The amplitude of the spectrum is given by 

| V(cjo) | 2 = | V 1 1 2 + | V 2 1 2 + 2V 1 V 2 coscd (t 2 -t 1 ) (3) 

2 

Examination of Eq . (3) shows that |v(cd)| varies sinusoidally 

2 2 

between the values (V + V ) and (V - V ) • Shaw and Gurnett [1971] 

I ^ 1 z 

used this model with V = and t - t^ = to explain some banding 

of whistlers observed on Injun 5. Their observations were made in the 
range 0 - 500 Hz and the spacings between bands typically ranged from 
10 to 30 Hz. These spacings correspond to a on the order of 50 

milliseconds. Uman [1969] reports maximum and minimum values (used in 
the sense that most values fall between these limits) for At of 100 
and 3 milliseconds, respectively with representative values of ~ 40 
milliseconds. The latter value fits Shaw and Gurnett's data and model 
very well, but the lower value is not sufficient to explain whistler "c" 
in Figure 3 which has a spacing of — 1300 Hz. A more serious objection 
is raised if one notes (as in Figure 3) that the spacing between nulls 
in the same banded whistler is often non-uniform. Furthermore, even 
among whistlers with uniform spacing (e.g., Figure 1), the nulls are 
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usually either much narrower than the maxima or vice-versa. However this 
model, based on Eq. (1), predicts a simple sinusoidal variation of 
amplitude with frequency rather than the spectral shapes noted above. 

For a discussion of the effects of substituting a more realistic model 
for the time behavior of a lightning stroke into Eq. (1) see the second 
paragraph in section on different paths within the ionosphere below. 

Lightning strokes from different locations . We can escape 
entirely the first objection raised above by assuming that the two 
strokes are not coming from the same lightning channel, but rather occur 
in different locations separated by a distance large enough to produce 
the observed range of frequency spacings. However, this model has precisely 
the same difficulty as the preceding one in terms of reproducing the 
observed spectral shapes . 

Different paths in the waveguide . A very similar tentative 
explanation is illustrated in Figure 11. There, one of the waves radiated 
by the lightning stroke travels directly to the entry point into the 
ionosphere while the other first reflects off the earth and hence travels 
a longer path. However, a Af ~ 1 kHz implies a travel time differential 
of ~ 1 millisecond corresponding to a path difference of ~ 300 kilometers. 
It is clear that the path differential maximizes when the entry point is 
directly above the source. Since the entry point is ~ 80 - 100 km above 
the earth's surface, not one, but several hops will be required. However, 
absorption and reflection losses would then make the multi-hop wave 
amplitude too small to cancel the direct wave to the extent shown in 
Figure 10. In addition this model is clearly no better able to reproduce 
the observed spectral shapes than the two preceding ones. 
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FIGURE 11. TENTATIVE EXPLANATION OF BANDING BY MULTIPLE PATH PROPAGATION IN THE ATMOSPHERE. 

The first path is a direct ray from the originating lightning source through 
the atmosphere to an entry point into the ionosphere. The second ray reflects off the 
earth’s surface while traveling between the same two end points. Higher order paths 
involving reflections off the ionosphere as well as the earth are also possible. For 
clarity the drawing is not to scale. 



Different paths within the ionosphere. The last possibility 


along the previous lines is that the nulls arise due to the superposition 
of two waves that have traveled different paths within the ionosphere. 
However, let us note (as, for example, in Figure 1) that the nulls may 
continue to be observed at a constant frequency for times on the order 
of a second. The diffuseness of these whistlers is attributed to multiple 
entry points into the ionosphere thereby implying multiple paths to the 
satellite, which is traveling at 7 km/sec . The phase-path integrals vary 
significantly along the paths, and it is unreasonable to expect that the 
same relative phase relationships will be maintined over a 7 km interval. 
Of course, this argument also applies to the two preceding models since 
At (and thus Af) will depend on the location of the entry point into 

the ionosphere in these two models. 

These explanations cannot be salvaged by assuming more realistic 
models for the time behavior of the lightning stroke in Eq. (1). Uman 
[1969] reports that the time duration (t g ) of a stroke is ~ 20 - 60 psec. 
Thus variations in a lightning pulse take place on a time scale < t g . It 
may then be shown by arguments similar to those on page 18 that these 
variations in the time domain produce spectral components which vary less 
rapidly than cos 2nf t . That is, t g would have to be one to two orders 
of magnitude larger than the accepted values to produce minima spaced as 
closely as those in the banded whistlers. 

3 . Banding in the Source Spectrum 

However, if such large times, t^, did actually exist, then we 
are saying in effect that the nulls arise from holes within the lightning 
source spectrum. If, indeed the sources were involved, then the same 
patterns should appear with similar regularity in the ground whistler 
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data. In fact, banding in the ground data appears, but only very rarely, 
and few examples have been found to date. 

4 . Gyroresonance Interactions 

The first few examples of banded whistlers noted by the author 
all had very uniform spacings . This immediately led to speculation that 
some type of gyroresonance interaction was involved. The existence of 
such interactions in the ionosphere was reported by Gurnett and Mosier 
[1969]. They found attenuation bands in auroral zone noise at integer 
multiples of the local proton gyrof requency . However, in Figure 1 the 
nulls occur at odd half multiples of 810 Hz (±30 Hz) rather than at 
integer multiples of the local proton gyrof requency of 719 Hz. Also in 
Figure 3 the local proton gyrof requencies range from 693 Hz at whistler 
"A" up to 728 Hz at whistler "e". Spacings observed in the latter whistler 
range from 430 Hz to 1000 Hz and little correlation seems to exist from 
one banded whistler to the next. For example, at "b" the spacing ranges 
from 300 Hz to 500 Hz while 18 seconds later at "c" a relatively uniform 
spacing of ~ 1300 Hz occurs. Therefore theories of ion gyroresonance 
interactions are not adequate to explain these observations. 

5 . Specular Reflection or Transmission 

The last model to be discussed in this report involves specular 
transmission or reflection of an incident plane wave from a slab of 
enhanced ionization aligned parallel to the earth T s magnetic field 
and extending in the E-W direction. Without loss of generality, B q 
is taken to be along the x-axis. The x-z plane is taken to be the earth's 
meridian plane, and the incident wave normal k is assumed to be in this 
plane. This is a common assumption, and as Thorne [1969] has pointed out, 
the medium has a tendency to focus k back onto the meridian plane if 
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the initial direction of k is not in this plane. Figure 12 shows the 
geometry . 

*■# 

Define: = incident electric field strength, 

E = reflected electric field strength, 

R 

= transmitted electric field strength, 

E = electric field of the upgoing wave in the slab 
+ 

E = electric field of the downoing wave in the slab 

= wave normal vector in region i 

k is related to the wave frequency od and the 
refractive index \i in the region by 

| k | = ^ | p. | where c is the velocity of light. 

E = electric field component in k direction 
n 

2^ = electric field component in X k(-y) direction 

E = electric field component in E X k direction 
n 

L = thickness of the slab. 

The boundary conditions are that the components of electric field 
parallel to a dielectric boundary are continuous across that boundary 
(see for example, Ramo et al., [1965]). Thus the conditions at the lower 
surface (assumed to be z = 0) are: 


E 

E 


XI 

YI 


+ E. 


XR 


= E. 


X+ 


+ E. 


YR 


= E. 


Y+ 


+ E 
4* E 


X- 

Y- 


The boundary conditions at the surface z = L are: 

ik L -ik L 

e xt = E E x + +E Z E x- 


ik L -xk L 

zB zB _ 

E = e E_. + e E 

YT Y+ Y- 


(4) 

(5) 

( 6 ) 

(7) 
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REGION A 



FIGURE 12, THE GEOMETRY OF THE WAVE-SLAB INTERACTION. 

In the upper figure, region A is the ambient 

plasma. Region B is a slab of thickness L whose density 

is enhanced by a factor 6 over that of the ambient plasma. 

The earth 1 s magnetic field B Q is in the x direction, 

parallel to the slab. The wave normal vector k lies in 

the plane formed by x and the slab normal direction, z. 

It makes an angle 8 with B q . The lower figure shows the 

three principal electric field axis of the wave, E^ 

(parallel to k), E t (in the -y direction), and E 

(in the E, x E direction), 
x n 
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Following Smith and Brice [1964] 


p = ±1,0 


2 _ , v r 

~ 1 2 pY + 1 

t' r r 

2 -2 

W = [i 
P P 

"e - 5 ( " 2 1 + W -l> 

2W 2 = 2W 2 + (W 2 -W 2 ) sin 2 © ± [(W 2 -W 2 ) 2 sin 4 © + (W 2 -W 2 ) 2 cos 2 ©] 4 ' 2 (8) 

~ ~ “ e o i-i 


o e 


where 


|a = refractive index 

W = l/p 

**■> 

9 = angle between k and 

k = k sinG evaluated in the slab 
zB 

£ = sum over all types of charged particles in the plasma 

r 

2 / 2 

X = f /f 
r or' 

Y = f„ /f 
r Hr 


f = wave frequency 
u Hr 


B q 

th ^ , 1 o 1 r 

f u = gyrof requency of r type particle = 2n m 


th 


f = plasma frequency of r type particle = 
or 


r 1/2 

N rV 

£ o M r 


1__ 

2 * 


c = dielectric constant of free space = 
o 

q 9 m , N are respectively the (signed) charge, mass, and 
r r r number density of the r th type particle. 

(rationalized MKS units used throughout) 

Note that 

EE E 

— = — cose - sine = Q 

EE E 

y y y 


Following Smith's class notes [1968]: 


009 90 9 9 222. COS0 

E w :E n :E = W 2 (W 2 V): - orV) (IT-iT) tan9: iW (W^W^) ~ r ~ 


(9) 
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Noting that E = - E we have 
y 1 

_(W 2 -w 2 )(W 2 -W 2 ) tan0 cos9 - W 2 (W 2 -W 2 ) sin9 
e o e 

Q — — 

2,2 2 . cos0 

-i W — — 


tan9 (W 2 -W 2 )W 2 

- 2i — ^ 5"“ (10) 

(w z -w ) w 


There is of course a different value of Q for each region and Q changes 
sign with 9 . Equations (4) and (6) transform to 


Q (E -E ) = Q (E -E ) 
YI YR u v Y+ Y- 


ik L -ik L 

E Y + ' e V>’ 


where the Q's are computed for the positive values of 9 


Using Eqs . (5) and (7) to eliminate E^ + and E^ 


-ik L 

E e Z (1 + y) - E (1 - y) = E (1 + y) 


+lk zB L 

E yt e (1 - Y ) - E yr (1 + y) = E yI (1 - v) 


where 


y = q a /q b 


Now eliminating E yields: 

YR 


def _ /_ 

T = E v /E. 


Vt/ y t 2 * 

2y cos k L - i(l + y ) sin k L 


Assuming no losses in the plasma (k and y are real) then the fraction 

Z 15 

of the wave energy transmitted through the slab (i.e., the transmission 
coefficient) is given by 


2 2 2 2 2 
4y cos kL+ (1+y ) sin (k L) 


±T 

2 22 . 2 

4y + (1“Y ) Sin (k zB L) 
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Under this assumption the reflection coefficient 
satisfies 


2 



1 



(15) 


Equation (14) requires knowledge of the value of the refractive 

index fj, and W (= [j. ) • In order to obtain this information Eq • (8) 

B B B 

must be combined with Snell's law: 


COS0 A = COS9 B 

which may also be written 

c°s 2 0 = a w 2 (16) 

B B 

where : 

a = cos 2 9 a /w| (17) 

Using Eq. (16) to eliminate 0 2 from Eq. (8) we obtain the following 

2 

quadratic expression in 

(l-aw 2 + aw 2 )w^ + (aw 2 w 2 -aw 2 w 2 -w 2 -w 2 )w 2 + w 2 w 2 = o (18) 

0 O B 1 —X COO tJ o c u 

where all quantities (except a) are evaluated inside the slab. Equation 

(18) (a particular form of Booker's quartic) is solved for the fast mode 

(i.e., choose the '+' sign in the solution) and then 0 B is determined 

with the aid of Eqs . (16) and (17). 0 , the slab thickness (L), and 

the quantities which determine W_^, W and (ion mass, f^, and 

f ) are all parameters. They are not necessarily independent. For 
H 

example if one specifies an altitude within the latitude range of 

occurrence of the banded whistlers this not only tends to specify the 

ions, f , and f„, but it may also, depending on one's model for 
P H 

slab formation, determine L. 
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To determine the transmission coefficients as a function of frequency 
a FORTRAN program (see Appendix) was written for use on Stanford Univer- 
sity’s IBM 360/67 computer. Several assumed models of slab thickness and 
altitude location were tried. These were combined with three different 
types of models for the variation of the incident wave normal angle, 9, 
with frequency. 

Before discussing these models let us note some of the properties 
of Eq. (14). This equation is deceptively simple and masks some rather 
complex physical and mathematical behavior. Calculations indicate that 
Y varies much more slowly than sin k L as f changes . To a good 
approximation the minimum values of Eq. (14) with respect to f are 
given by 

T | 2 = ~^~2~2 <19> 

min (1 + y ) 

This quantity is plotted in Figure 13 as a function of the incident wave 

normal angle 6 . As expected it goes to zero at zero wave normal angle 

since the wave is then traveling parallel to the slab. As the wave normal 

angle increases, more and more energy is directed into the slab. However, 

at the same time the impedance mismatch between the two media starts to 

grow as the difference in refractive indices increases. Thus, it is not 

unreasonable to expect that at some point the latter effect might counter 

balance the former and It I . would reach a peak. 

* min 

The slab also has two other interesting properties. However, we 
must first distinguish between the three cases illustrated in Figure 14. 
Figure 14a sketches the closed wave normal surfaces for regions A and B 
characteristic of frequencies below the lower hybrid resonance (LHR) . 

The approximate expression for the LHR frequency was given by Brice and 
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30 40 50 60 70 80 90 


20 30 40 50 


FIGURE 13. VARIATION OF THE MINIMUM VALUE OF THE TRANSMISSION COEFFICIENT , | T | - n (AS 
DEFINED IN EQ . (19) IN THE TEXT), WITH THE INCIDENT WAVE NORMAL ANG1JE, Qj • 
The graphs are parameterized by the wave frequency f in kiloHertz. The following 
parameters for the ambient medium (A) were used in the calculations: plasma frequency 

f =3.2 MHz, electron gyrof requency f = 1.3 MHz and ionic composition 100% 0 + , 
leading to a lower hybrid resonance at 7703 kHz. The left graph corresponds to an 
enhancement factor 6=5 and the right one to 6=2. 




and Smith [1965] as: 


LHR 


,2 ~ 2 

1 f O 1 H 

M 2 2 

f o + f H 


1/2 


( 20 ) 


where M is the ion to electron mass ratio of the single ion (0 ) 
assumed in this model. Since we are dealing with a slab of enhanced 
ionization it is obvious that there will be a range of frequencies that 
lie above the LHR frequency for region A and below that for region B. 

In this range the refractive index surfaces will appear as shown in 
Figure 14b. Above this range both surfaces open up as shown in Figure 
14c. It is interesting to note in the latter two cases that the two 
surfaces cross . 

Figure 13 also illustrates that for each frequency below the LHR 
frequency there is an angle at which perfect transmission occurs regard- 
less of slab thickness. (It happens when y = 1 — Eq . (19) — i.e., 
when the ratio E E is the same in both media.) This phenomenon is 
still unexplained. It does not appear to be the usual Brewster angle 
[Ramo et al., 1965] effect since in the present case the electric field 
of the wave is not polarized in the plane of incidence. 

Finally, for frequencies corresponding to Figure 14c there exists 
a range of angles for which no real solution of Snell's law exists in 
the slab. Due to the finite thickness of the slab some tunneling may 
be expected to occur. However, assuming that the slab is several wave- 
lengths thick, it becomes very strongly reflecting. The behavior of the 
10 kHz lines in Figure 13 aptly illustrates this phenomenon. 


B. COMPUTATIONAL RESULTS 

The simplest and perhaps most instructive model for the banded 
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FIGURE 14. SKETCH OF REFRACTIVE INDEX SURFACES [i AND \± B FOR 
MEDIA A AND B RESPECTIVELY. a) At frequencies 
below the lower hybrid resonance (LHR) both surfaces are closed, 
b) At frequencies above the LHR in region A but below that in 
region B the surface in region A opens and intersects that of 
region B. c) At frequencies above the LHR in both regions both 
surfaces are open and intersect. As indicated by the dotted line 
there is a range of incident wave normal angles for which no 
solution in region B exists. (For clarity these drawings are not 

to scale.) 
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whistler phenomenon is to assume that waves of all frequencies are 
incident upon the slab at the same wave normal angle. First, note that 
expression (14) for |t| 2 maximizes whenever |sin(k zB L)| = 0 that is 

f |i B (f,0g) sin0 B = — m = 1,2,... (21) 


For sufficiently small values of 9 , we may express the refractive 

B 

index in the quasi-longitudinal (QL) approximation (see Helliwell [1965]) 

Then, p 2 ~ x/y cos0 . Thus the lefthand side of Eq . (21) will vary as 

2 

/f implying that f is proportional to m . This means that the 
transmission maxima (or reflection minima) will get progressively 
farther apart as f increases. At the other extreme, let us assume 
that the constant angle 0 D is so large that it is near the resonance 
cone angle corresponding to a frequency of 10 kHz. Then as f -♦ 10 kHz, 
jj. oo . Thus, the successive frequency solutions of Eq . (21) get pro- 
gressively closer together. It then becomes intuitively plausible that 
somewhere between these two extremes there will exist angles that will 
produce a relatively uniform spacing. Figure 15 shows the variation of 
reflection coefficient with frequency for an incident wave normal angle 
of 80°, and it confirms the above expectation. 


Recognizing the highly idealized nature of this model, let us 
examine some of its consequences. First, and perhaps most important, 
the maxima in Figure 15 are much broader than the minima. This differs 
radically from the spectral amplitude variations predicted by some of 
the models discussed earlier, and it is in agreement with the spectral 
shapes observed on many banded whistlers such as the one in Figure 1. 

However, this model cannot account for spectra such as that in 
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whistler 'A' of Figure 3. As indicated above, small incident wave 
normal angles cause a progressive increase in the spacing of the maxima, 
large angles produce a decrease, and a small range of intermediate 
angles produce a relatively uniform spacing. There are no angles at 
which highly irregular spacings occur. 



FIGURE 15. REFLECTION SPECTRUM. The variation of the 

reflection coefficient with frequency is shown. 
The ambient plasma has one ion (0 + ) a plasma frequency of 
3.2 MHz, and an electron gyrof requency of 1.3 MHz. The 
slab has a density five times ambient and a thickness of 
0.35 km. Waves at all frequencies are assumed to be 
incident at an angle of 10° from the normal to the slab. 
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Finally, at high wave normal angles, Figure 13 indicates that 
large density enhancement factors (say > 5) would be required to produce 
the observed attenuations of — 14 db in the transmitted spectra. If 
we are dealing with reflected spectra, then the situation is slightly 
different. While large attenuation is easily obtained (the reflected 
signal amplitude goes to zero when Eq . (21) is satisfied), the problem 
now is how to reflect enough of the wave to produce a reasonably strong 
signal without requiring large enhancement factors. While irregularities 
of this size are known to exist in the auroral regions [Lund et al., 

1967] no such variations have been reported at the latitudes where banded 
whistlers occur. In fact, most observers (see the review article by 
Calvert and Warnock [1969]) report variations normally on the order of 
a few percent at midlatitudes. 

In summary the assumption that waves of all frequencies are incident 
on the slab at the same angle requires first, large wave normal angles, 
and second, high density enhancement factors in order to account for a 
whistler such as that in Figure 1. However, this assumption is not 
necessary, and is probably not even tenable. At the very least different 
frequencies are refracted by different amounts upon entering the ionosphere 
with the spread in angles proportional to the square root of the frequency. 
However, for a plasma frequency of 3 MHz and electron gyrof requency of 
1 MHz this will only produce a maximum angular spread on the order of 
1-2 degrees in the frequency range of 1-10 kHz. Horizontal gradients 
may act to increase this spread. 

The objections raised in the preceding three paragraphs can be 
considerably reduced or eliminated by using much smaller incident wave 
normal angles that are however strongly frequency dependent. This latter 
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assumption appears to be quite reasonable. For example, ray tracing 
calculations by Scarabucci [1969] have indicated that the various 
frequency components of a 1-hop whistler must have a large variation 
in wave normal angles in order that they all be able to focus on a given 
satellite position. Rough calculations based on his Figures 45 and 49 
indicate a rise on the order of some 30° in the wave normal angles with 
respect to the earth's magnetic field for a 1-hop whistler at 120 km 
altitude in the frequency range of 2-10 kHz with sharp drop-offs above 
and below this range. The variation could be modeled either by a saw- 
tooth or a sine curve. It must be borne in mind that Scarabucci 's 
results were highly dependent on his assumption of a model for electron 
density variation. Nevertheless, it appears reasonable to conclude that 
the wave normal angle of a wave component is highly dependent on its 
frequency . 

A variety of models for the variation of the incident angle 0 
with f were tried. Figure 16 illustrates a highly irregular trans- 
mission spectrum. Note the qualitative resemblance to whistler 'A' in 
Figure 3. Figure 17 shows a reflection spectrum for the same slab, but 
with a different model for G^f). Note the exceptionally uniform 
spacing of the nulls in this spectrum. There is, of course, nothing 
unique about this model. 

The procedure for reproducing a particular banded whistler's 

spectrum is relatively straightforward once the plasma parameters (i.e., 

slab location) f and f have been selected. Using the QL approxi- 
p H 

mation, calculation of the wavelength in the slab (using a "reasonable" 
value for the enhancement factor, 6) allows calculation of a fairly 
good estimate for the thickness L which will produce a spacing that 
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0 2 4 6 8 10 

f, kHz 


FIGURE 16. COMPUTED TRANSMISSION COEFFICIENTS. The 

ambient plasma has one ion (0 + ), plasma frequency 
of 3.2 MHz, and electron gyrof requency of 1.3 MHz. The slab 
has 2.0 km thickness and a density twice that of the ambient. 
The variation of incident wave normal angle 0^ with fre- 
quency shown above produces an irregular spacing of the maxima. 
Note how the peaks first have a spacing of ~ 1 kHz, then 
~ 1.8 kHz, then ~ 1 kHz again. 


will on the average approximate that observed in the whistler . A computer 
run (taking ~ 2 seconds) assuming that the incident angle 9^ increases 
linearly by 30° over the 1-10 kHz interval is made. From the results L 
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may be adjusted if the initial average spacing is too far off. Then 
adjustments are made to 9j(f) bearing in mind that increasing 0j(f) 
in a given frequency interval acts to decrease the spacing between nulls 
in that range and vice-versa. Figure 17 illustrates a calculated reflection 
spectrum in which a very uniform spacing between minima was obtained by 
suitable variation of incident wave normal angle 9^. with frequency. 

If amplitude information on the whistler is available to be matched, small 
changes in 6 will not affect the calculated spectral shape if the 
product L 6 is held fixed. 
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FIGURE 17. COMPUTED REFLECTION COEFFICIENTS. With the 
same parameters as in Figure 16, but with a 
different variation of incident wave normal angles with 
frequency, a very uniform spacing of ~ 1600 Hz is produced. 
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The exact details of how the wave arrives at the slab and what 


the actual dependence of 0 on frequency is require more detailed ray 
tracing calculations. Two plausible models are schematically illustrated 
in Figure 18. 
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MODEL I 


MODEL n 



FIGURE 18. TWO WAVE-SLAB INTERACTION MODELS— MODEL I APPEARS TO BE USABLE ON ALL 1-HOP 
BANDED WHISTLERS. Model II is required for 0 + banded whistlers. It can 
only produce monotonically varying spacing on 1-hop whistlers. With small enhancement 
factors and small wave normal angles it can easily reproduce a uniform spacing. Trans- 
mission is, of course, also possible for both models. Interaction from the outside of 
the field line is considered less likely since waves tend to head outward toward higher 
L values. 



IV. DISCUSSION 


In the previous chapter several possible mechanisms to explain the 
banded whistler phenomenon were examined. Only the wave-slab interaction 
model provided positive results in terms of explaining the observed 
spectral shapes of the banded whistlers. 

This model may also find application in explaining the phenomena 
shown in Figures 2a and 2b. For example, if the slab were not sharp 
sided, the ratio of maxima to minima in the spectrum can be expected to 
reduce. We might well expect to obtain spectra resembling that in 
Figure 2b. 

If this mechanism is indeed the correct one, then it may provide 
us with a valuable tool for studying some types of ionospheric irregu- 
larities. More data should be examined to determine how the banded 
whistlers (and, by implication, the irregularities) vary with magnetic 
storm activity (preliminary results indicate that little, if any, 
correlation exists) and local time. 

A particularly good correlation already appears to exist with 
the data on HF and MF ducting. In a rocket-borne HF sounding experiment 
Knecht et al., [1961] reported irregularities with thicknesses of about 
1 to 10 km. Loftus et al., [1966] using the Explorer 20 topside-sounder 
found thicknesses ranging from 1 to 40 km. Average thicknesses of 
0.6 km were observed by Muldrew [1963] using Alouette 1 data. These are 
consistent with the thicknesses used in the preceding chapter. In 
addition, Ramasastry and Walsh [1969] report that ducting is seen up to 
regions near the plasmapause (while Muldrew [1967] reports that the 
ducts cut off at L = 1.3 to 1.4) which are consistent with the data in 
Chapter 2 . 
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Another avenue for further research is to improve the rigor of the 
calculations in the preceding chapter. For completeness, the non-whistler 
mode should also be included in Eqs . (4) through (7) (providing four new 
variables) and the conditions for continuity of magnetic fields across 
di-electric boundaries should be used (thus providing the necessary two 
new equations at each boundary). For the frequencies of interest in 
this report, these waves are evanescent, decaying exponentially with 
distance away from the slab boundaries where they may be excited. To 
say that these waves are identically zero or have zero "skin depth" 

[Ramo et al . , 1965] is equivalent to saying that the medium behaves like 
a perfect conductor to this mode and is thus capable of supporting surface 
currents. In the presence of surface currents the tangential B fields 
are not continuous across the surface. Thus the choice of continuity 
of E rather than B fields in Eqs. (4) through (7) seems consistent 
with the assumption that the non-whistler node waves could be neglected. 
Hence, we can see that this assumption has profound mathematical impli- 
cations. However, from the physical point of view, it is not expected 
that the greatly increased mathematical complexity which will result from 
relaxing this assumption will significantly change the results because 
these extra waves are nonpropagating . 

It should also be noted that Eq . (8) ignored collisions. This is 

2 

justified by the fact that the collision rate, which is of order 10 at 
satellite altitudes, is small compared with the whistler frequencies. 

In summary, the mechanism of spectral reflection or transmission 
from a field aligned slab of enhanced ionization appears to be best able 
to explain the observed spectral shapes of the banded whistlers. The 
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model assumes that the slab's density is much larger (at least by a 
factor of two) tha that of the ambient plasma. Since such irregularities 
have not yet been reported in the literature as having been seen at the 
mid-latitudes where the banded whistlers occur, this assumption must 
stand as a prediction of the model. 
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APPENDIX 


: PrOG'-’AM SLAB 

E. r: A | M ASS 

D I MR MS I PM Fl5?r»» t T(50C) 
c OM MOV MA S S , F P , r H , W P , WL » W F , WO 
P!=B. 14159265^6 

RAO=PI/lRn. 

MAS$=2 0?76. 

F ° = ? ? 0 0 0 0 0 . 

FH = 1 360000 • 
r>tLTA= 1. 5 
ANGL F= 16 . 

THKNS=4. 

F{ 1 ) = 1 OOO. 
r>r-TN T 510 

51 i' FORMAT ( * FREQUENCY ANGLE R "LFCY TOM COEF !C TENTS *, / *<H7*> 

HO 500 J =1,201 

CALCULATE THE SQUARE OF THF PHASE VELOCITY IN THE AMBIENT PLASMA 
PHI =P( I ,F { J) , ANGLE , PI, RAD) 

CALL FORMOL) ( F ( J ) 1 
A= ( WF - W H ) * S T M ( P H I) **? 

W1=WF-. 5*A+SQRT( A*A + ( ( WR- Wl ) *COS( PHI ) ) **? ) *. 5 
IF(Wl) 502, 50.?, ^50 

45 0 P0LAR1 =?.*TAN( PH! ) * { HE - W 1 ) * ( WO/ W 1 ) / ( W R -W L I 
C * CALCULATE THE SQUARE OF THE PHASE VELOCITY IN THE SLAB 
FP=FP*SQRT ( OELT A) 

CALL FOP MOD (F ( J) 1 
FP = F°/ S31K0ELTA ) 

D=C OS l PMI >**2/Wl 
A = 1 . -D*wF+D*WO 
B = 0*( WR*WL-WE*WO) -wo-wf 
BP=B#P-4. * A*WO*WE 
I F { B 5 ) 502,470, 4 7 Q 

W?= .5* ( -R+SOR T ( PP ) ) / A 
IF(W2*D-1.) 490, ^°D, 5 0 2 
49 0 PHI 2=ARC0S( SQF T( W2*0H 
c**»* CALCULATE THE REFLECTION CPcFFICISNT 

X=( ( TAN( PHI 2 ) *2 . * ( WF-W2 ) * (W 0/W2 ) / ( WR-WL M / POL &R1 ) **2 
WA VE= 1 . 5F 5* SORT ( W ? ) /( -( J ) *P T *SI N( PHI 2 ) ) 

A = S I N ( THKNS/W AVF ) **2 

T ( J ) = 1 . -4 . * X / ( 4 • * X* ( 1 • - A ) 4- A* ( 1 • + X ) * * 2 ) 

fcC I «= ( J - T ) 600, 600, 501 

6(0 WA VFL=2.*WAVE*P I \ 

56 1 CONTINUE 
E = PHI/P.AO 

PRINT 5 1 ? , F ( J ) , E , T ( J ) 
p FORMAT ( « * , F 10 . 0 , F 6. 2 , 6 !. 6. 6 ) 

500 F(JH ) =F ( J 1+50. 

502 PRINT 60B , WAVE L 

6C ° FORM AT (• WAVELENGTH AT 1 *H7 TS E16.6) 

STOP 

E f iD 



43 


SEL-71-054 



S UB E ni J T I ■| c F 0 0 w P n ( F ) 

C**** TMI S routine co v putes thf rquar c of the phase vel 3 CI t i fs of t h f 

r>:** * LIMITING MODES I V J A ONE ION! PLASMA. WR AND WL REPRESENT THE 
r**** VALUES AT ZERO PFGRFES OF THE FAST AMD SLOW MODES RESPECTIVELY. 
£ m##* SIMILARLY, WF AMD WO APE THE RESPECTIVE VALUES EVALUATED AT 
(•*-*** NINc T Y DEGREES. 

P r A L MASS 

CO’-' MON MASS,FP,pH»WR , W L , W r » W *. 3 

X= ( F n / F ) 

V-FH/f 

WE=l./I l.-X*( MAS S«-l.)/m. -Yl*c MASS* Y)M 
a'L = 1 . / { 1 .-X*( MASS + I . )/ ( ( 1 . +Y ) *( MASS-Y ) ) ) 

wr=i. z( l.-x-x/MASSi 

WF = ( 1 . +X * ( M ASS + I • ) * ( Y* Y-M A c S) / ( ( Y* Y - 1 • ) *( Y^Y-VASS**?) ) ) * WF * WL 

RETURN 

END 


Ft ) M C T I ON P( I,G, H, PI, PAD) 

THIS POUT INF COMPUTES THF VALUF OF THF: INCIDENT AN5LE»P* AS A 
- **** FUNCTION of THf FREQUENCY .0 , AND THF INCIDENT ANSI E AT 1 KHZ* H. 
IF (G- 4000. ) 1,1,? 

] P=( H-( G-1000. > /?DQ. I *PAD 

RETURN 

? IF ( G-6500 • ) 3,3,'+ 

3 P=( H-15.M G-AOOG. ) * 7. /?500 • ) *RAD 

R c T I J R N 

a p=(H- R . +• ( G-6 5 C 0* • ) * 1 2 • / A 5 00 • ) *R A D 
RETURN 
c\ir) 
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